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Abstract 


Shortly after its emergence in southern China in 2002/2003, Severe Acute Respiratory 
Syndrome coronavirus (SARS-CoV) was confirmed to be the cause of SARS. Subsequently, 
SARS-related CoVs (SARSr-CoVs) were found in palm civets from live animal markets in 
Guangdong and in various horseshoe bat species, which were believed to be the ultimate 
reservoir of SARSr-CoV. Till November 2018, 313 SARSr-CoV genomes have been 
sequenced, including 274 from human, 18 from civets and 47 from bats [mostly from 
Chinese horseshoe bats ( Rhinolophus sinicus), n=30; and greater horseshoe bats 
(Rhinolophus ferrumequinum), n=9]. Hie human SARS-CoVs and civet SARSr-CoVs were 
collected in 2003/2004, while bat SARSr-CoVs were continuously isolated in the past 13 
years even after the cessation of the SARS epidemic. SARSr-CoVs belong to the subgenus 
Scirbecovirus (previously lineage B) of genus Betacoronavirus and occupy a unique 
phylogenetic position. Overall, it is observed that the SARSr-CoV genomes from bats in 
Yunnan province of China possess the highest nucleotide identity to those from civets. It is 
evident from both multiple alignment and phylogenetic analyses that some genes of a 
particular SARSr-CoV from bats may possess higher while other genes possess much lower 
nucleotide identity to the corresponding genes of SARSr-CoV from human/civets, resulting 
in the shift of phylogenetic position in different phylogenetic trees. Our current model on the 
origin of SARS is that the human SARS-CoV that caused the epidemic in 2002/2003 was 
probably a result of multiple recombination events from a number of SARSr-CoV ancestors 
in different horseshoe bat species. 
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1. Introduction 


Shortly after its emergence in southern China in 2002/2003, Severe Acute Respiratory 
Syndrome coronavims (SARS-CoV) was confirmed to be the cause of SARS (Peiris et al., 
2003). There has been a total of 8,096 laboratory-confirmed cases of SARS, leading to 774 
mortalities in 11 countries (World Health Organization, 2004). Subsequently, SARS-related 
CoV (SARSr-CoV) was found in palm civets from live animal markets in Guangdong 
province (Guan et al, 2003). However, the initial hypothesis that civets may act as animal 
reservoir of SARS-CoV was soon overturned by several observations. Firstly, SARSr-CoV 
was only detected in civets from the market, but not those in the wild (Kan et al, 2005; Tu et 
al, 2004). Secondly, the high ratio of nonsynonymous to synonymous mutation rates (Ka/Ks 
ratios) of the Spike (S), open reading frame (ORF) 3a and non-structural protein (nsp) genes 
in civet SARSr-CoVs collected in both the 2003 and the minor 2004 outbreaks suggested that 
the vims was undergoing rapid evolutionary gene adaptation in civets (Song et aL, 2005). 
Thirdly, compared with SARS-CoV collected from human during the 2003 epidemic, 
functional changes have been observed in the S protein of civet SARSr-CoV and the SARS- 
CoV isolated from the 2004 minor outbreak. Hie latter showed less efficient use of the 
human angiotensin converting enzyme 2 (ACE2) receptor (Li et al, 2005b) and demonstrated 
resistance to antibody inhibition (Yang et al., 2005). Finally, while significant levels of 
antibody to SARS-CoV were detected in 80% of the civets from one animal market in 
Guangzhou, low seroprevalence rates in civets from various civet farms in China suggested 
that civets were largely brought to the animal market infection-free (Guan et al, 2003; Tu et 
al, 2004). The SARSr-CoV was likely contracted during later mixing and trading. 

In view of these observations, we carried out a molecular surveillance study in various 
mammals in Hong Kong to hunt for the ultimate source of the vims (Lau et aL, 2005). 
Among the 127 bats (including 8 bat species), 20 monkeys and 60 rodents surveyed, SARSr- 




CoV was only detected in 39% of the fecal samples of 59 Chinese horseshoe bats 


(Rhinolophus sinicus, Rs) (Lau et aL, 2005). Western blot analysis showed that antibodies 
against the nucleocapsid (N) protein of bat SARSr-CoV was present in 67% of the serum 
samples of Chinese horseshoe bats, while 8% of the serum of Chinese horseshoe bats were 
tested positive for human SARS-CoV-neutralizing antibody with titer > 1:20 (Lau et aL, 
2005). Shortly afterward, another independent group also reported the detection of SARSr- 
CoV in Chinese horseshoe bats, greater horseshoe bats (Rhinolophus ferrumequinum, Rl), 
and big-eared horseshoe bats (Rhinolophus macrotis, Rm) in Hubei and Guangxi provinces of 
China (Li et aL, 2005a). In the past few years, SARSr-CoVs have been isolated from a 
variety of horseshoe bats in Yunnan province of China by several groups (Ge et aL, 2013; He 
et aL, 2014; Hu et aL, 2017; Lau et aL, 2015). 

Three hundred and thirteen SARSr-CoV genomes have been sequenced from 2003 to 2018. 
These include 274 genomes from human, 18 from civets and 47 from bats. Hie human 
SARS-CoV and civet SARSr-CoV were collected in 2003/2004, while bat SARSr-CoVs 
were continuously detected even after the cessation of the SARS epidemic. In this article, we 
review our current understanding of the molecular epidemiology, evolution and, phylogeny 
of SARSr-CoVs based on analysis of these 313 genomes. 




2. The SARS-CoV genome 


Hie genome size of the SARS-CoV varies from 29.0 kb to 30.2 kb. Its genome structure 
follows the characteristic gene order of other known CoVs: the 5’ two thirds of the genome 
comprises ORFlab encoding replicase polyproteins, while the 3’ one third consists of genes 
encoding structural proteins including S, envelope (E), membrane (M), and N proteins (Fig. 
1). Both the 5’ and 3’ ends of the SARS-CoV genome contain short untranslated regions. The 
translational product of ORFlab is cleaved by proteases encoded by SARS-CoV itself into 16 
nsps, which include major enzymes such as papain-like protease(s) (PF pro ), chymotrypsin- 
like protease (3CF pro ), RNA-dependent RNA polymerase (RdRp) and helicase (Hel) (Fig. 1). 
In contrast to the genome of viruses belonging to lineage A Betacoronavirus (recently 
renamed as subgenera Embecovirus), haemagglutinin-esterase gene is absent from the 
genome of SARS-CoV. In addition, SARS-CoV contains 6-7 accessory proteins, encoded by 
ORF3a, ORF3b, ORF6, ORF7a, ORF7b and ORF8 [or ORF8a and ORF8b as a result of a 29- 
nucleotide (nt) deletion]. This is unique to lineage B Betacoronavirus , a sub genus recently 
renamed Sarbecovirus and contains all SARSr-CoVs (Fig. 2). 

Studies in the past 15 years have partly revealed the biochemical functions of these accessory 
proteins (Fiu et al., 2014). Protein 3a triggers apoptosis and induces the production of 
proinflammatory cytokines such as RANTES (Regulated on Activation, Normal T cell 
Expressed and Secreted; also known as C-C motif chemokine ligand 5, CCF5) and CXCF8 
(C-X-C motif chemokine ligand 8). Protein 3b inhibits type I interferon and also induces 
apoptosis. Protein 6 inhibits interferon signaling and stimulates DNA synthesis. Protein 7a 
activates NF-kB (nuclear factor kappa B) and MAPK8 (mitogen-activated protein kinase 8) 
for CXCF8 and RANTES production. Hie function of protein 7b is not well characterized 
yet. ORF8 is present in all SARSr-CoV genomes in bats and civets, as well as in SARS-CoVs 
isolated from human during the early phase of the epidemic. Protein 8 activates the ATF6 




(activating transcription factor 6) branch of unfolded protein response. In the genomes of 
SARS-CoVs isolated from human during the late phase of the epidemic, there was a signature 
29-nt deletion in ORF8, splitting it into two separate ORFs 8a and 8b. Protein 8a includes 
caspase-dependent apoptosis whereas protein 8b modulates cellular DNA synthesis. 




3. Unique phylogenetic position of SARS-CoV 


Before the SARS epidemic, there were just around 10 CoVs with complete genome 
sequences available. These Co Vs were classified into three groups: Group 1, Group 2 and 
Group 3. In 2011, the Coronavirus Study Group of the International Committee for 
Taxonomy of Viruses renamed these three groups into three genera: Alphacoronavirus, 
Betacoronavirus, and Gammacoronavirus (de Groot et al, 2011). When SARS-CoV was first 
discovered in 2003, phylogenetic analysis of the SARS-CoV genome showed that it occupied 
a unique position in Betacoronavirus , which was subsequently placed into the subgenus 
Sarbecovirus. The traditional betaCoVs (e.g. mouse hepatitis virus, human CoV OC43, 
bovine CoV) were classified as Embecovirus (Fig. 2). After the SARS epidemic, an 
unprecedented number of novel CoVs were discovered (Lau et al, 2016; Lau et al., 2007; 
Lau et aL, 2012a; Lau et aL, 2012b; Lau et aL, 2014; Woo et aL, 2005; Woo et aL, 2014a; 
Woo et al, 2014b). This led to the description of lineage C Betacoronavirus , which 
comprises important members such as Tylonycteris bat coronavirus HKU4, Pipistrellus bat 
coronavirus HKU5, Hypsugo bat coronavirus HKU25 and Middle East Respiratory 
Syndrome CoV (Lau et al., 2013; Lau et aL, 2018b; Woo et al., 2007; Woo et al., 2006a), and 
lineage D Betacoronavirus (Lau et al., 2010b) as well as a novel genus Deltacoronavirus 
(Lau et al, 2018a; Woo et al, 2012; Woo et al, 2017) (Fig. 2). Lineage C and lineage D 
Betacoronavirus were now renamed as sub genera Merbecovirus and Nobecovirus. 




4. Molecular epidemiology and evolution of SARS-CoV 

4.1. Circulation of SARSr-CoV in horseshoe bats in 2004 to 2018 

Since its first discovery in Chinese horseshoe bats in 2004 (Lau et aL, 2005), SARSr-CoVs 
have been continuously found in various horseshoe bat species in the last 13 years (Drexler et 
aL, 2010; Ge et aL, 2013; He et aL, 2014; Hu et aL, 2017; Lau et al., 2005; Lau et aL, 2010a; 
Li et aL, 2005a; Tang et al., 2006; Wu et al., 2016; Yang et aL, 2013; Yang et aL, 2015; Zeng 
et aL, 2016). This is in contrast to the case for civets and human, where SARSr-CoVs were 
only found in 2003/2004, and never reported afterward. For the 47 bat SARSr-CoV genomes, 
30 are from Chinese horseshoe bats, 9 from greater horseshoe bats, 2 from big-eared 
horseshoe bat, 2 from least horseshoe bat ( Rhinolophus pusillus, Rp), and 1 each from 
intermediate horseshoe bat ( Rhinolophus affinis, Ra), Blasius’s horseshoe bat ( Rhinolophus 
blasii, Rb), Stoliczka’s Asian trident bat ( Aselliscus stoliczkanus. As) in the neighboring 
family Hipposideridae and wrinkled-lipped free-tailed bat ( Chaerephon plicata, Cp) in the 
genetically more distant family Molossidae. SARSr-CoVs have also been detected in 
countries other than China, including Thailand, Italy, Luxembourg, Bulgaria, Slovenia, 
Hungary, Japan, Kenya, etc. However, only partial sequences were available for these 
isolates. Nevertheless, the immediate progenitor of SARS-CoV has not been pinpointed. 

4.2. Geographical gradient of SARSr-CoV 

Among the 45 SARSr-CoV genomes from bats in China, 11 of them were from Hong Kong, 
2 were from Guangdong, 2 from Guangxi, 5 from Hubei, 20 from Yunnan and one each from 
Shaanxi, Shanxi, Jilin, Guizhou, Hebei, respectively. Overall, it is observed that the SARSr- 
CoV genomes from bats in Yunnan possess relatively higher, whereas those from 
Guangdong, Jilin, Shanxi and Hebei possess relatively lower nt identity to those from civets 


(Fig. 3A). This is an interesting phenomenon since the SARS epidemic that emerged in late 




2002 was first noticed in Guangdong province as an outbreak of acute community-acquired 
atypical pneumonia syndrome. Severe cases were later retrospectively traced back to five 
cities around Guangzhou, with the index case reported in Foshan, a city 24 km away from 
Guangzhou. The second case was a chef from Heyuan who worked in a restaurant in 
Shenzhen who contacted with wild game-food animals regularly. This mismatch between 
clinical events and apparent gradient of nt identity could either be due to a missing link 
during evolutionary adaptation among the SARSr-CoVs in different provinces or simply as a 
result of sampling error. 

4.3. Recombination and evolution 

The high frequency of homologous RNA recombination is one of the major factors 
contributing to a plastic genome underpinning the evolutionary force in CoVs. This has 
resulted in different genotypes or even different CoVs adapted to new hosts (Herrewegh et 
aL, 1998; Lau et al., 2011; Terada et aL, 2014; Woo et al., 2006b). As for SARSr-CoVs, it is 
evident from both multiple alignments and phylogenetic analyses that some genes of a 
particular SARSr-CoV from bats may possess higher while other genes possess much lower 
nt identity to the corresponding genes of SARSr-CoV genomes from human/civets, resulting 
in shifting of phylogenetic position in different phylogenetic frees. This phenomenon is 
frequently observed in SARSr-CoVs and likely explains the generation of novel SARSr- 
CoVs that could jump from bat to civet and subsequently to human. 

4.4. Sprotein of SARSr-BatCoVs 

Trimers of S protein form spikes on the surface of CoVs particles. It comprises two 
functionally distinct subunits - SI and S2 domains which are involved in receptor binding 
and fusion respectively. I . ike other class I viral fusion proteins, the S protein undergoes a 




series of events including receptor recognition, proteolytic cleavage to shed the SI subunit 
and conformational changes in S2 that ultimately lead to fusion of the viral and host 
membranes. It has been well established that human SARS-CoV utilizes ACE2 as a 
functional receptor. The receptor-binding motifs in the C-terminal domain of SI are 
implicated in receptor recognition. Substitutions within the SI receptor-binding domain 
(RBD) confers adaptability to new or orthologous entry receptors, thus altering the viral 
tropism (Hulswil et aL, 2016). Undoubtedly, the ability to bind human ACE2 is an 
indispensable step in establishing cross-species transmission. 

As we have previously hypothesized, the five amino acid (a.a.) deletion, twelve a.a. deletion 
as well as a.a. substitutions at 5 critical residues for binding serve as determining factors for 
the S protein-ACE2 interaction (Table 1). A critical pre-requisite for reliable ACE2 
utilization is the absence of 5 a.a. and 12 a.a. deletions and, preferably, the presence of at 
least two out of five human-adapted residues. Based on these analyses, SARSr-Rs-BatCoV 
WIV1 is one of the strains most advantageously conformed genotypically for ACE2 
utilization. It has been shown to be able to directly infect well-differentiated primary human 
airway epithelial cell cultures (Menachery et aL, 2016). In addition, neutralization assays 
using convalescent sera from SARS-CoV infected patients showed robust neutralization 
against tissue culture infectious dose 50ofWIVl (Ge et al., 2013). 

SARSr-Rs-BatCoV RsSHC014, a strain discovered in Chinese horseshoe bats in 2012, 
contains two of the five a.a. residues in civet strain civet007 but none of the five human- 
adapted residues. It retains both the 5 a.a. and 12 a.a. deletion sites. This genotype is shared 
by SARSr-Rs-BatCoV Rs4231 and Rs4084. Recombinant mouse-adapted SARS-CoV 
expressing the S protein of RsSHC014 was still able to utilize ACE2 for viral entry, causing 
cytopathic changes in Vero cells and weight loss in mice model, despite the apparent failure 


in pseudovirus infectivily assay (Menachery et aL, 2015). However, significantly slower viral 




replication rate was observed, suggesting that deletions in RBDs were more critical in 
receptor recognition while the presence of ACE2-adapted critical residues modulated entry 
efficiency. In contrary to the WIV1 strain, neither neutralizing human monoclonal antibody 
nor existing double-inactivated whole SARS-CoV vaccine provided protective effect against 
infection caused by SHC014-harboring recombinant vims strain, indicating key difference in 
a.a. sequences determining antigenicity (Menachery et aL, 2015). 

A phylogenetic tree was constructed from the RBDs of SARSr-CoVs (Fig. 4). Surprisingly, 
SARSr-Ra-BatCoV LYRall, a strain discovered in intermediate horseshoe bat in Baoshan, 
Yunnan (He et aL, 2014), around 375 km from Kunming, sits closest to the civet/human 
SARSr-CoVs despite its lower overall nt identity. This is likely due to the presence of 7 nt in 
civet/human SARS-CoV that is found exclusively in SARSr-Ra-BatCoV LYRall. Such a 
unique S gene in relation to civet/human SARS-CoV would therefore challenge the previous 
hypothesis that the origin of RBD was solely from Chinese horseshoe bats in Kunming. 
Nevertheless, the S proteins of the majority of bat SARSr-CoVs found in Chinese horseshoe 
bats in various part of China including Yunnan resemble the genotype of SARSr-Rs-BatCoV 
HKU3 but not the human ACE2-utilizing genotype of SARSr-Rs-BatCoV WIVl. This raises 
question on how SARSr-BatCoVs manage to have two distinct genotypes of S infecting the 
same host (Rs). There is a possibility that two cellular receptors for SARSr-CoVs are present 
in Chinese horseshoe bats. 

Interestingly, two novel SARSr-Rp-BatCoV strains ZXC21 and ZC45 isolated from least 
horseshoe bat in Zhoushan, Zhejiang province in eastern China were shown to be able to 
infect suckling rats, causing inflammation in the brain tissue and histological changes in the 
lung and intestine despite failed viral isolation in VeroE6 cells. Phylogenetically, these two 
strains represented a separate clade, lying between the susceptible clade of SARSr-Ra- 


BatCoV LYRall and the non-susceptible clade of SARSr-Rf-BatCoV Rfl. However, the 




mechanism of such infectivity in suckling rats was not investigated further, especially in light 
of the two retained deletion sites in the S genes of ZXC21 and ZC45, which should have 
imposed inter-species barrier in terms of receptor specificity based on previous discussion 
(Hu etaL, 2018). 

4.5. SARSr-CoV in Chinese horseshoe bats and greater horseshoe bats 
Among all the horseshoe bats, SARSr-CoVs are most commonly found in Chinese horseshoe 
bats and greater horseshoe bats. Detailed analysis of their genomes revealed several 
intriguing phenomena. Firstly, ORFla and ORFlb in SARSr-Rf-BatCoV YNLF_31C 
isolated from greater horseshoe bats in Lufeng, Yunnan by our group has the highest nt 
identity to civet SARSr-CoV, especially at the regions nsp5, nsplO, nspl2 and nspl3 (Fig. 
3A). This raised the possibility that Lufeng, Yunnan could be a potential gene pool for the 
progenitor of ORF1 of SARS-CoV. Secondly, as previously illustrated, for most regions 
along the genome of SARSr-CoV except ORF8, SARSr-Rs-BatCoVs from Chinese 
horseshoe bats are predominantly closer to civet SARSr-CoV. However, the ORF8 of 
SARSr-Rf-BatCoV has higher a.a. identity to civet SARSr-CoV ORF8. In addition, V90 and 
1113 are the two Rf-specific residues present in the ORF8 of SARSr-Rf-BatCoV Rf4092, 
which could evolutionarily bridge the ORF8 of SARSr-Rf-BatCoV to SARSr-Rs-BatCoV 
that had very high identity to civet/human SARS-CoV (Fig. 5). SARS-Rf-BatCoV was thus 
hypothesized to be the origin of human SARS-CoV ORF8 gene through recombination 
events. 

In fact, at least three genotypes of ORF8 have been found to circulate in nature (Wu et al, 
2016). Type I is genetically closest to ORF8 of civet SARSr-CoV (96.2-98.1% nt identity). 
There are eleven such strains isolated from bats so far, with eight of them originating from 


Chinese horseshoe bats and three from greater horseshoe bats. Type II ORF8 has lower (70.8- 




82.8%) nt identity to ORF8 of civet SARSr-CoV. It is prevalent only in SARSr-Rf-BatCoV 
and demonstrates genetic stability regardless of geographical distribution with similar ORF8s 
detected in Korea, Shanxi, Hubei, and Hebei SARSr-Rs-BatCoV has not been found to 
process type II ORF8. Yet another strain, the SARSr-Rp-BatCoV F46 isolated from least 
horseshoe bats, processes an ORF8 with a phylogenetic position right between type I and 
type II ORF8s with 93.2% nt identity to civet SARSr-CoV, further complicating the analysis 
(Fig. 6). Hie rest of the ORF8 genes from sequenced SARSr-BatCoV strains belong to type 
IE. These type HI ORF8 genes possess 53.9-57.7% nt identity to the ORF8 of civet SARSr- 
CoV and have the greatest genetic distance observed over the whole genome of SARSr-CoV 
(Fig. 3B). They are detected across a wide range of bat species including R. sinicus, R. 
macrotis, R pusillus, R. affinis, A. stoliczkanus and C. plicata, except R. ferrumequinum in 
which type HI ORF8 has not been found. In short, the ORF8 of SARSr-Rf-BatCoV from 
greater horseshoe bats is phylogenetically closer to ORF8 of human/civet SARS-CoV. 

As mentioned above, ORF8 in the late-phase human epidemic was split into 8a and 8b due to 
a unique 29 nt deletion. A similar 5-nt deletion was observed in the type I ORF8 isolated 
from the strain SARSr-Rs-BatCoV Rs4084. However, since the majority of SARSr-Rs- 
BatCoVs’ ORF8 genes belong to the type III ORF8 genotype, it raises the speculation that 
this minority of type I ORF8 plays a non-essential role in the Chinese horseshoe bats host 
environment. 

In a recent study, recombinant SARS-CoV carrying three forms of ORF8 were constructed, 
including complete ORF8 as seen in the early phase of the 2003 human epidemic and civets, 
ORF8 that contained the characteristic 29 nt deletion, and a variant strain with ORF8 
replaced by a 5 nt spacer sequence. Replication assay in VeroFM cells demonstrated a 
phenotypic hierarchy that the infectious clone with full type I ORF8 exhibited significantly 


more efficient growth than other variants. The same hierarchy of viral replication was 




observed in human airway epithelial cell cultures and several “non-host” cell lines transduced 


to express human ACE2, including R. alcyone lung epithelial cell line, cotton rat airway 
epithelial cell line, and goat and sheep lung cell lines. This study suggests that the 29 nt 
deletion actually conferred a loss of viral fitness during the initial phase of human-to-human 
transmission, contradicting the previous belief that the truncated products ORF8a and ORF8b 
favored adaptation to human (Muth et al., 2018). 

4.6. Genomic overview of SARSr-BatCoVs 

Complete genome analysis has been carried out for all sequenced SARSr-CoV with respect to 
civet SARSr-CoV. The importance of civet SARSr-CoV is underlined by its close genetic 
relationship to human SARS-CoV along the entire genome and especially at the RBD region 
where nt identity approaches 98.8-99.8%. Moreover, the genetic diversity within the group of 
civet SARSr-CoV strains is minimal, further supporting that civets are only recently infected. 
Adaptation of SARSr-BatCoV to civet as intermediate amplification host is likely needed, 
which is further corroborated by the phylogenetic position of civet ACE2 in-between bat and 
human ACE2. Recently, direct bat-to-human transmission of SARSr-CoV was proposed by 
some groups in view of the fact that some SARSr-BatCoVs could utilize human ACE2 
receptors in vitro and that seropositive sera against SARSr-BatCoVs were detected in 
Yunnan residents (Menachery et aL, 2015; Wang et al., 2018), despite the fact that 
human/civet SARS-CoV has never been reported in Yunnan. Our findings provide substantial 
evidence against this hypothesis from a phyloepidemiologic point of view and demonstrated 
the significance of an intermediate host during the transmission cascade (Fig. 3). 

Referring to Fig. 3A, SARSr-Rs-BatCoV Rs4084 was one of the two strains that possess the 
highest genome identity to civet/human SARSr-CoV. It is so far the only SARSr-BatCoV 
that possessed at least 70% nt identity in all parts of the genome, including the four 




hypervariable regions. The region sharing lowest similarity to civet/human SARSr-CoV was 
the N-terminal domain (NTD) of the S protein which was likely non-critical As mentioned 
above, SARSr-Rs-BatCoV Rs4084 was predicted to be able to utilize human ACE2 as 
functional receptor since its RBD region is genetically similar to RsSHC014. SARSr-Rs- 
BatCoV Rs4084 has not been studied either as whole vims or in the form of recombinant 
virus in terms of accessory protein function, tissue infectivity assay and replication kinetics; 
leaving a potential field for future investigations. Animal experiment using civet model 
would offer great promise in the studying of such vims, along with other SARSr-BatCoVs 
that could utilize human ACE2. 




5. Concluding remarks 


Hie >300 genome sequences of SARS-CoV in human, civets and bats accumulated in the last 
15 years allowed us to understand the cause of such a highly fatal epidemic that affected us in 
a global scale in 2003. After 15 years of work, we are now much closer to fully understand 
the origin of SARS-CoV and its evolution, as genomes with fragments that contain gene 
sequences with higher and higher nt identities to the SARS-CoV found in human were 
observed. Our current model on the origin of SARS is that the human SARS-CoV that caused 
the epidemic in 2003 was probably a result of multiple recombination events from a number 
of SARS-CoV ancestors. Yunnan is the province in China with the largest diversity of 
horseshoe bat species. It is also the area shown to have a high variety of SARS-CoV 
ancestors. Further in-depth molecular epidemiology studies in this and other provinces in its 
proximity, including Guangxi and Guangdong, will hopefully give us an even clearer picture 
on the origin and early evolution of SARS-CoV. 
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Legends to figures 


Fig. 1. Genome organization of SARS-CoV. ORFlab with nspl-16 are colored in blue. 
Structural proteins including S, E, M and N are in pink. Accessory proteins were numbered 
and in yellow. 

Fig. 2. Hie diversity of Co Vs as demonstrated with a phylogenetic tree targeting the 303bp 
partial RdRp sequence (position 15293-15596 with respect to Human SARS-CoV TOR2). 
The neighbor-joining phylogenetic tree was constructed with maximum composite likelihood 
method by MEGA 7.0. The test of phytogeny was statistically supported by the bootstraps 
value calculated from 1,000 frees. Bat Co Vs were labeled with black triangles. Hie branches 
of Gammacoronavirus and Deltacoronavirus were compressed. 

Fig. 3. (A) Nucleotide identity and (B) Genetic Distance along the SARS-CoV genome for 
all available SARSr-CoVs with complete genomes. A comparison is made with reference to 
Civet SARSr-CoV SZ3. Strains are listed in (A) descending order from the top according the 
whole genome identity (B) ascending order from the top according the whole genome genetic 
distance at the last column. Red and blue boxes represent the highest and lowest end of the 
identity, respectively. White boxes represent deletions. The titles of SARSr-Rf-BatCoVs 
from greater horseshoe bats are highlighted in pink. Hie data were generated by Matrix 
Global Alignment Tool (MatGAT) (Campanella et aL, 2003). 

Fig. 4. Phylogenetic tree constructed based on the nt sequences of the RBD of the S protein 
of SARS-CoVs and SARSr-CoV. SARSr-Rf-BatCoVs are labeled with dots. Brackets in red, 
orange, yellow and blue represented the descending order of nt identity with respect to civet 
SARSr-CoV SZ16. Hie phylogenetic tree was constructed by Maximum L ikelihood method 
with T92+G as the substitution model by MEGA 6.0. Hie test of phytogeny was statistically 
supported by the bootstraps value calculated from 1,000 trees. 




Fig. 5. Multiple alignment of type I ORF8 of SARSr-Rs-BatCoV YN2013, GX2013, civet 
SARSr-CoV and SARSr-Rf-BatCoV strain Rf4092 against the type II ORF8 from other 
SARSr-Rf-BatCoVs. Host-specific residues were highlighted with red boxes. 

Fig. 6. Phylogenetic tree constructed based on the nt sequences of ORF8 of SARS-CoVs. 
SARSr-Rf-BatCoVs are labeled with dots. Brackets in red, orange, yellow and blue 
represented descending order of nt identity with respect to civet SARSr-CoV SZ16. The 
phylogenetic free was constructed by Maximum L ikelihood method with T92+G as the 
substitution model by MEGA 6.0. The test of phytogeny was statistically supported by the 
bootstraps value calculated from 1,000 trees. 




Table 1. Summary of the critical elements for ACE2 utilization present in the S protein of 
various SARSr-BatCoVs 



NTD 
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12 a.a. 
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SARSr-BatCoVs that could replicate in cell lines ■ 

SARSr-BatCoVs that were reported to have slower replication kinetics (Menachery et al. 2015) 
SARSr-BatCoVs that failed to replicate in cell lines 

# Critical a.a. residues on S protein determining interaction with ACE2 

* Only RBD sequences are available in Genbank 








Highlights 


• 313 SARSr-CoV genomes have been sequenced (274 from human, 18 civets and 47 
bats) 

• SARSr-CoV genomes of bats in Yunnan possess highest nt identity to those from 
civets 

• The origin of human SARS-CoV was probably a result of multiple recombination 
events 

• Recombination from a number of SARSr-CoV ancestors in different horseshoe bat 
species 
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SARSr-Rf-BatCoV 273/2005/R.ferrumequinum/Hubei/2005 

SARSr-Rs-BatCoV HKU3-12/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-6/R.sinicus/Hong Kong/2005 
SARSr-Rs-8atCoV HKU3-4/R.sinicus<Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-2/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-5/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-10/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-11/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-1/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-9/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-3/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-13/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-8/R.sinicus/Guangdong/2006 
SARSr-Rs-BatCoV HKU3-7/R.sinicus/Guangdong/2006 
SARSr-Rf-BatCoV SX2013/R.lerrumequinum/Shanxi/2013 

SARSr-Rf-BatCoV HeB2013/R.ferrumequinum/Hebei/2013 
SARSr-Rf-BatCoV JL2012/R.»emimequinum/Jilin/2012 

Rf-BatCoV 16BO133/R.fefrumequinum/Korea/2016 


SARSr-Rb-BatCoV BM48-31/R.blasii/Bulgaria/2008 


Figure 3A 

























































SARSr-CoV strains 

Civet SARSr-CoV SZ16/Shenzhen/2003 
Human SARS-CoV GZ02/Guangzhou/Feb2003 
Human SARS-CoV TOR2/Toronto/Mar2003 
Human SARS-CoV BJOI/Beijing 
Civet SARSr-CoV civet010/Guangdong/2004 
Civet SARSr-CoV B039/Guangdong/2004 
Civet SARSr-CoV GD69/Guangdong/May2003 
Civet SARSr-CoV civet007/Guangdong/2004 
Civet SARSr-CoV A022/Guangdong/2004 



Civet SARSr-CoV civet020/Guangdong/2004 
Civet SARSr-CoV PC4-227/Guangdong/2004 
Human SARSr-CoV GZ0401/Guangzhou/Dec2003 
Civet SARSr-CoV PC4-136/Guangdong/2004 
Civet SARSr-CoV PC4-13/Guangdong/2004 
Civet SARSr-CoV A001 
Human SARSr-CoV GZ0402 
SARSr-Rs-BatCoV WIV16/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4874/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4231/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs9401/R.sinicus/Yunnan/2015 
SARSr-Rs-BatCoV Rs7327/R.sinicus/Yunnan/2014 
SARSr-Rs-BatCoV Rs3367/R.sinicus/Yunnan/2012 


SARSr-Rs-BatCoV WIV1 

SARSr-Rs-BatCoV RsSHC014/R.sinicus/Yunnan/2011 
SARSr-Rs-BatCoV Rs4084/R.sinicus/Yunnan/2012 
SARSr-As-BatCoV As6526/Aselliscus stoliczkanus/Yunnan/2014 
SARSr-Rs-BatCoV Rs4237/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4255/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4247/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV YN2013/R.sinicus/Y unnan/2013 
SARSr-Rs-BatCoV Rs4081/R.sinicus/Yunnan/2012 

SARSr-Rs-BatCoV Rs672/R.sini, u-, <,uizh. ,uVOdr, 

SARSr-Rf-BatCoV YNLF~3 

SARSr-Rs-BatCoV Rp3/R.sinicus/Guangxi/2004 
SARSr-Rs-BatCoV GX2013/R.sinicus/Guangxi/2013 
UNVERIFIED: SARSr-Rp-BatCoV F46/R.pusillusA'unnan/2012 
S ARSr-Ra-BatCoV LYRal 1 /R.affinis/Yunnan/2011 
SARSr-Cp-BatCoV Yunnan2011/Chaerephon plicataA'unnan/2011 
SARSr-Rs-BatCoV HuB2013/R.sinicus/Hubei/2013 
SARSr-Rp-BatCoV Shaanxi2011/R.pusillus/Shaanxi/2011 
SARSr-Rm-BatCoV 279/200SIR.macrotis/Hubei/2005 
SARSr-Rm-BatCoV Rml/R.macrotis/Hubei/2004 
SARSr-Rs-BatCoV HKU3-2/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-4/R.sinicus/Hona Kona/2005 
SARSr-Rs-BatCoV HKU3-1/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-3/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-6/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-5/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-12/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-9/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-10/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-11/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-13/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-7/R.sinicus/Guangdong/2006 
SARSr-Rs-BatCoV HKU3-8/R.sinicus/Guangdong/2006 


SARSr-Rf-BatCoV SX2013/R.femimequinum/Shanxi/2013 
SARSr-Rf-8atCoV JL2012/R.fenumequinum/Jilin/2012 
SARSr-Rf-BatCoV 16BO133/R.ferrumequinum/Korea/2016 

SARSr-Rb-BatC oV BM48-31 /R.blasii/Bu lgaria/2008 




Figure 3B 

















Civet SARSr-CoV PC4-13/Guangdong/2004 
Civet SARSr-CoV A022/Guangdong/2004 
Civet SARSr-CoV civet010/Guangdong/2004 
Civet SARSr-CoV civet007/Guangdong/2004 
Civet SARSr-CoVcivet020/Guangdong/2004 
Civet SARSr-CoV A001 
Civet SARSr-CoV PC4-136/Guangdong/2004 
Civet SARSr-CoV B039/Guang dong/2004 
Civet SARSr-CoV PC4-227/Guangdong/2004 
Human SARSr-CoV GZ0401/Guangzhou/Dec2003 
Human SARSr-CoV GZ0402 
Civet SARSr-CoV SZ3/Shenzhen/2003 
Civet SARSr-CoV SZ16/Shenzhen/2003 
Human SARS-CoV GZ02/Guangzhou/Feb2003 
Human SARS-CoV BJ01/Beijing 
Civet SARSr-CoVGD69/Guangdong/May2003 
Human SARS-CoV TOR2/Toronto/Mar2003 
— SARSr-Ra-BatCoV LYRall/R.affinis/Yunnan/2011 
SARSr-Rs-BatCoV WlVI6/R. sinicus/Y unnan/2013 
SARSr-Rs-BatCoV Rs4874/R. sinicus/Y unnan/2013 
SARSr-Rs-BatCoV WIV1 

1 SARSr-Rs-BatCoV Rs3367/R.sinicus/Yunnan/2012 
SARSr-Rs-BatCoV Rs7327/R.sinicus/Yunnan/2014 
SARSr-Rs-BatCoV Rs9401/R. sinicus/Y unnan/2015 
SARSr- Rs-BatCoV Rs4231/R.sinic us/Y unna n/2013 
SARSr-Rs-BatCoV Rs4084/R. sinicus/Y unnan/2012 
SARSr-Rs-BatCoV RsSHC014/R.sinicus/Yunnan/2011 

-SARSr-Rb-BatCoV BM48-31/R.blasii/Bulgaria/2008 


SARSr-Rs-BatCoV HuB2013/R.sinicus/Hubei/2013 

SARSr-Rs-BatCoV HKU3-3/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-2/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-1/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-4/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-5/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-6/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-9/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-10/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-13/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-11/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-12/R.sinicus/Hong Kong/2007 
, SARSr-Rs-BatCoV HKU3-8/R.sinicus/Guangdong/2006 
ipo I SARSr-Rs-BatCoV HKU3-7/R,sinicus/Guangdong/2006 
SARSr-Rs-BatCoV GX2013/R.sinicus/Guangxi/2013 
— SARSr-Rs-BatCoV Rs4247/R. sinicus/Y unnan/2013 
— SARSr-Rs-BatCoV Rp3/R.sinicus/Guangxi/2004 
93 . SARSr-Rm-BatCoV Rml/R.macrotis/Hubei/2004 
—ftSARSr-Rm-BatCoV 279/2005/R.macrotis/Hubei/2005 

' -SARSr-Rp-BatCoV Shaanxi2011/R.pusillus/Shaanxi/2011 

SARSr-Cp-BatCoV Yunnan2011/Chaerephon plicata/Yunnan/2011 
Hr— SARSr-As-BatCoV As6526/Aselliscus stoliczkanus/Yunnan/2014 
SARSr-Rs-BatCoV Rs4255/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4237/R.sinicus/Yunnan/2013 

■ SARSr-Rs-BatCoV Rs672/R.sinicus/Guizhou/2006 

■ SARSr-Rs-BatCoV Rs4081/R.sinicus/Yunnan/2012 
- SARSr-Rs-BatCoV YN2013/R.sinicus/Yunnan/2013 
100 r# SARSr-Rf-BatCoV Rf4092/R.ferrumequinum/Yunnan/2012 

UNVERIFIED: SARSr-Rp-BatCoV F46/R.pusillus/China/2012 
SARSr-Rf-BatCoV Rfl/R.ferrumequinum/Hubei/2004 
i SARSr-Rf-BatCoV 273/2005/R.ferrumequinum/Hubei/2005 
SARSr-Rf-BatCoV SX2013/R.ferrumequinum/Shanxi/2013 
SARSr-Rf-BatCoV HeB2013/R.ferrumequinum/Hebei/2013 
I SARSr-Rf-BatCoV YNLF 31C/R.ferrumequinum/Yunnan/2013 
I SARSr-Rf-BatCoV YNLF 34C/R.ferrumequinum/Yunnan/2013 
SARSr-Rf-BatCoV 16BO133/R.ferrumequinum/Korea/2016 
90 1 V SARSr-Rf-BatCoV JL2012/R.ferrumequinum/Jilin/2012 


L_ 81 . 


99 


99 


00 1 1 

42 


o.i 


Figure 4 






































20 


30 


40 


50 


60 


Civet SARS-CoV SZ3 
SARSr-Rs-BatCoV YN2013 
SARSr-Rs-BatCoV GX2013 


10 

- I -|- | - | -|- | - | -|- | - | -|- | 

MKLLIVLTCISLCSCIRTWQRCASNKPHVLEDPCPTGYQPEWNIRYNTRGNTYSTAWLC 

. H . 


SARSr-Rf-BatCoV 

Rfl 

. .F. 

. .M. . . 

. . .C. 


.T. 

.Y. . 

. .N. . . 


.NS R. . 

SARSr-Rf-BatCoV 

273/2005 

. .F. , 

. .M. . . 

. . .C. 


.T, 

.Y. . 

. .N. . . 


.NS R. . 

SARSr-Rf-BatCoV 

YNLF 34C 



. . .C. 

. V. 

.T, 

.Y. . 

. .N. . . 

. . R. 

.N. R. . 

SARSr-Rf-BatCoV 

YNLF 31C 



. . .C. 

. V. 

.T. 

.Y. . 

. .N. . . 

. . . R . 

.N. R. . 

SARSr-Rf-BatCoV 

SX2013 


. .M. . . 

. . .C. 

. V. 

.T, 

.Y. . 

. .N. . . 


.N. R. . 

SARSr-Rf-BatCoV 

HeB2013 


..M. . . 

. . .C. 

. V. 

.T, 

.Y. . 

. .N. . . 


.N. R. . 

SARSr-Rp-BatCoV 

F46 







. .N. . . 


. K.G . 


Civet SARS-CoV SZ3 
SARSr-Rs-BatCoV YN2013 
SARSr-Rs-BatCoV GX2013 
SARSr-Rf-BatCoV Rf4092 
SARSr-Rf-BatCoV Rfl 
SARSr-Rf-BatCoV 273/2005 
SARSr-Rf-BatCoV YNLF_34C 
SARSr-Rf-BatCoV YNLF_31C 
SARSr-Rf-BatCoV SX2013 
SARSr-Rf-BatCoV HeB2013 
SARSr-Rp-BatCoV F46 


70 80 9C 

- I - | - | - | - | . . . | 

ALGKVLPFHRWHTMVQTCTPNVTINCQD] 7! 


100 110 

I-I- I -I. 

!G LI RCWYLHEGHQTAAFRI 


FYQ.P 

FYQ. S 
F.K.S 
F.K. S 
FYRCS 

FYR. S 


120 


I 


VtWLT-KRT 
. N- . 

. N- . 

. N- . 

1 . D . FF . 

I. D . FF . 

I . D . FF . 
i.D.FF. 

I. D . FF . 
i . D . FF . 


Figure 5 



























































Human SARSr-CoV GZ0401/Guangzhou/Dec2003 
Human SARSr-CoV GZ0402 
Civet SARSr-CoV civet020/Guangdong/2004 
Civet SARSr-CoV A001 
Civet SARSr-CoV A022/Guangdong/2004 
Civet SARSr-CoV civet007/Guangdong/2004 
Civet SARSr-CoV civet010/Guangdong/2004 
Civet SARSr-CoV B039/Guangdong/2004 
Civet SARSr-CoV PC4-227/Guangdong/2004 
Civet SARSr-CoV PC4-13/Guangdong/2004 
Civet SARSr-CoV PC4-136/Guangdong/2004 
Civet SARSr-CoV SZ16/Shenzhen/2003 
Civet SARSr-CoV GD69/Guangdong/May2003 
Human SARS-CoVTOR2/Toronto/Mar2003 
Human SARS-CoV BJ01/Beijing 
Civet SARSr-CoV SZ3/Shenzhen/2003 
Human SARS-CoV GZ02/Guangzhou/Feb2003 

— 6 SARSr-Rf-BatCoV 16BO133/R.ferrumequinum/Korea/2016 
I SARSr-Rf-BatCoV SX2013/R.ferrumequinum/Shanxi/2013 
SARSr-Rf-BatCoV HeB2013/R.ferrumequinum/Hebei/2013 
(SARSr-Rf-BatCoV Rfl/R.ferrumequinum/Hubei/2004 
i SARSr-Rf-BatCoV 273/2005/R.ferrumequinum/Hubei/2005 
SARSr-Rf-BatCoV YNLF 31C/R.ferrumequinum/Yunnan/2013 
93 ! # SARSr-Rf-BatCoV YNLF 34C/R.ferrumequinum/Yunnan/2013 

— UNVERIFIED: SARSr-Rp-BatCoV F46/R.pusillus/Yunnan/2012 
SARSr-Rs-BatCoV Rs4084/R.sinicus/Yunnan/2012 
# SARSr-Rf-BatCoV Rf4092/R.ferrumequinum/Yunnan/2012 

- SARSr-Rs-BatCoV GX2013/R.sinicus/Guangxi/2013 

SARSr-Rs-BatCoV YN2013/R.sinicus/Yunnan/2013 
99 j SARSr-Rm-BatCoV 279/2005/R.macrotis/Hubei/2005 
83 J1 SARSr-Rm-BatCoV Rml/R.macrotis/Hubei/2004 
7A 1-SARSr-Rs-BatCoV HuB2013/R.sinicus/Hubei/2013 

— SARSr-Rs-BatCoV Rp3/R.sinicus/Guangxi/2004 
— SARSr-Rp-BatCoV Shaanxi2011/R.pusillus/Shaanxi/2011 
99 j SARSr-Ra-BatCoV LYRal 1/R.affinis/Yunnan/2011 

SARSr-Cp-BatCoV Yunnan2011/Chaerephon plicata/Yunnan/2011 
SARSr-Rs-BatCoV Rs7327/R.sinicus/Yu nnan/2014 
SARSr-Rs-BatCoV WIV1 

SARSr-Rs-BatCoV Rs3367/R.sinicus/Yunnan/2012 
SARSr-Rs-BatCoV Rs4237/R.sinicus/Yunnan/2013 
SARSr-As-BatCoV As6526/Aselliscus stoliczkanus/Yunnan/2014 
SARSr-Rs-BatCoV Rs4255/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4081/R.sinicus/Yunnan/2012 
SARSr-Rs-BatCoV RsSHC014/R.sinicus/Yunnan/2011 
- SARSr-Rs-BatCoV Rs672/R.sinicus/Guizhou/2006 
SARSr-Rs-BatCoV Rs4247/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs9401/R.sinicus/Yunnan/2015 
SARSr-Rs-BatCoV Rs4874/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV WIV16/R.sinicus/Yunnan/2013 
SARSr-Rs-BatCoV Rs4231/R.sinicus/Yunnan/2013 

-SARSr-Rs-BatCoV HKU3-8/R.sinicus/Guangdong/2006 

SARSr-Rs-BatCoV HKU3-7/R.sinicus/Guangdong/2006 
SARSr-Rs-BatCoV HKU3-13/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-6/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-4/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-12/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-11/R.sinicus/Hong Kong/2007 
SARSr-Rs-BatCoV HKU3-10/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-9/R.sinicus/Hong Kong/2006 
SARSr-Rs-BatCoV HKU3-5/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-1/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-3/R.sinicus/Hong Kong/2005 
SARSr-Rs-BatCoV HKU3-2/R.sinicus/Hong Kong/2005 


Figure 6 

























